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We investigate the three flavor Nambu-Jona Lasinio model of neutral quark matter at zero tem- 
perature and finite density, keeping into account the scalar, the pseudoscalar and the Kobayashi- 
Maskawa-'t Hooft interactions as well as the repulsive vector plus axial- vector interaction terms 
(vector extended NJL, VENJL in the following). We focus on the effect of the vector interaction on 
the chiral restoration at finite density in neutral matter. We also study the evolution of the charged 
pseudoscalar meson energies as a function of the quark chemical potential. 



I. INTRODUCTION 



Quantum Chromodynamics (QCD) is accepted nowadays as the theory of strong interactions. The phase diagram 
of strongly interacting matter in the T — fi plane, where T , \x denote respectively the temperature and the baryon 
chemical potential, is one of the most intriguing research topics in modern theoretical physics. Investigations on the 
various phases of the QCD phase diagram enlighten the physics of high temperature/small density, as well as small 
O ■ temperature/large density, matter: the former can be created in our laboratories by means of heavy ions collisions; 
the latter may be realized in the cores of the compact stellar objects (white dwarfs and neutron stars). The equation 
of state of strongly interacting matter is all what one needs to make theoretical predictions in the whole T — fx 
regime. Unfortunately, QCD is analytically treatable only in the pcrturbative regime. Moreover, the most interesting 
phenomena arc not of a pcrturbative nature. For this reason, the QCD equation of state is only known in a small 
slice of the phase diagram. 

The most concrete knowledge about strong interactions at /i = comes from Lattice QCD (LQCD). However, at 
fi > LQCD calculations with three colors suffer from the sign problem. In order to circumvent this problem several 
approaches have been suggested: expansion in [ijT P, 0, reweighting tecniques Q and analytical continuation from 
the imaginary chemical potential axes 0, [![. However none of them has not yet been of a practical use at high 



' chemical potential and small temperature region. Therefore, to make theoretical investigations on the QCD phase 



diagram, and to compute the equation of state in a wider range of chemical potentials where perturbative QCD does 
■ not work, some effective models are needed. Among them the Nambu-Jona Lasinio model 0, 0, H, [1] is probably 
Q\ | the most popular one: it shares the global symmetries of QCD and is simpler to handle than QCD itself. The main 
characteristic of the QCD vacuum, that is spontaneous breaking of the chiral symmetry, is described in the NJL 
model in a clear way. Moreover its minimal extension to the Polyakov-NJL model [l(J, [Tl| reproduces LQCD results 
at [i, = as well as at small chemical potential [l2| and at imaginary chemical potential [r|. As a consequence NJL 
' model is a promising tool to make calculations of the QCD phase diagram. 

In this work, we investigate the three flavor Nambu-Jona Lasinio model of neutral and /3-cquilibrated quark matter 
at zero temperature and finite density, keeping into account the scalar, the pseudoscalar and the Kobayashi-Maskawa- 
't Hooft interactions [HI, EH as well as the repulsive vector plus axial- vector interaction terms [HI, [13, EH Ell ■ The 
introduction of the vector interaction and thus of the vector excitations is interesting for several reasons: firstly, it 
is well known that vector interactions play a dominant role in determining the properties of matter at intermediate 
densities. Experimental progresses on the measurements of the in-medium properties of vector mesons can be found 
in Refs. [26|, [23|. In these conditions of baryon densities the main degrees of freedom are nucleons and the lighter 
pseudoscalar mesons (it and K), interacting among them by the exchange of vector mesons (mainly p and cj). One 
may wish to derive meson-baryon interactions starting from a microscopic interacting theory of the constituent 
quarks. Such a derivation can be obtained by means of the well known bosonization and hadronization of the NJL 
lagrangian [20l . l2l| . Keeping in mind the relevance of vector meson modes on the phenomenology of nuclear matter at 
intermediate densities, the vector interaction (which excites vector and pscudovector mesons) must be included from 
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the very beginning in the hadronization of the NJL lagrangian. Secondly, vector meson exchange might be responsible 
for kaon condensation at high density. With respect to this phenomenon, it is of a certain interest to compare the 
scenario offered by the NJL model with that obtained within the Hidden Local Symmetry framework, the latter being 
consistent with kaon condensation [2^, HH, HH, HH ■ Thirdly, it has been suggested that quark hadron continuity can 
be realized by means of the spectral function continuity of the vector mesons [28j . 

These are only a part of the reasons that lead us to consider the role of the vector interaction in quark matter, 
with particular focus on neutral and (3— equilibrated systems which should be realized in the core of neutron stars, if 
it exists. 

We study in this work the extended NJL model with vector interaction, with particular emphasis put on the 
neutral and equilibrated quark matter. The role of the vector interaction in NJL model has been studied several 
times 0, [l?! O, d, [H, H3, [3l| , but not yet in the regime of electrical neutrality and (3 equilibrium. The study of 
neutral, as well as /3-equilibrated, matter is interesting for astrophysical applications: matter inside a compact star 
has to be electrically neutral; moreover, since weak processes have a small characteristic time compared to the lifetime 
of a compact star, matter inside it should be equilibrated with respect to weak-interactions also. We can anticipate 
one of the results of our study, namely the phase diagram we compute does not differ qualitatively from that found 
in the aforementioned references. However, this does not diminish the relevance of our work. As a matter of fact, if 
one wishes to study the ground state of the model, having in mind physical applications as for example the structure 
of the compact stellar objects, then neutrality is a necessary requirement that must be settled in, and not a mere 
academic problem. 

The vector coupling Gy in the vacuum can be determined by the fit of the vector meson spectrum. However it is 



not clear if Gy in the nuclear medium has the same value it has in the vacuum [31[ . Experimental data on the medium 
modifications of the w vector meson mass [27| for densities up to the saturation nuclear density can be reproduced 
within a NJL model with a fixed value of Gy but with additive multiquark interactions (32l . |33| ; the effect of these 
interactions can be rephrased simply as density dependent redefinitions of the coupling constants [32l . l33l . l34l . |35| . As 
a consequence it is not wrong to think to a model in which only four quark interactions are included but where the 
coupling constants run with density. Instead of computing Gy at each value of the chemical potential in this work we 
treat it as a free parameter, and investigate on its influence on the restoration of the approximate chiral symmetry 
at finite density. NJL studies with Gy = predict a first order chiral restoration at p of the order of the constituent 
quark masses. When Gy is switched on, if its magnitude is larger than a critical value, then the transition becomes a 
smooth crossover [l9| . We reproduce this scenario with neutral matter at equilibrium. Also, a strong enough vector 
interaction disfavors the existence of the critical end point of the phase diagram [3l[ . 

We do not include color superconductivity in this study. This is done for simplicity. Nevertheless our results for 
the electron chemical potential and for the densities of the various species bear some implication of the effects of the 
repulsive vector interaction in the color superconductive phases. A quantitative study on the role of vector (as well as 
multiquark) interaction on the two flavor color superconductor has been performed in Refs. (30l.[35j. but an analogous 
study in realistic three flavor quark matter is still missing. 

Another point that we consider in some detail is the spectrum of the pseudoscalar excitations. This subject has 
been considered, within the NJL model but in the non-neutral case, in Ref. [36j . We firstly derive the loop expansion 
of the meson action. It includes scalar, pseudoscalar, vector and pseudovector mesons. However we focus here only on 
the pseudoscalar modes, therefore we specialize our equations to this case. We leave a more complete study to a next 
work. This choice is motivated by the comparison we would be tempted to make with nuclear matter models, which 
predict a kaon condensed phase for densities in the range (2.5 — 5)po H3]- We anticipate another result, namely the 
absence of strangeness condensation even in presence of a vector interaction, although the vector interaction lowers 
the K~ in-medium energy. This result was anticipated in Ref. [3(j| . However, the authors of this paper do not make 
any attempt to compute the chemical potential felt by the kaons, which is a relevant ingredient to exclude definitely 
meson condensation in the ground state, since they are not interested to the neutral phases. Moreover, they made 
their analysis on the basis of the low density approximation within the lowest linear level in p/po- Even though at 
this order some model independent prediction is possible, one might be rather interested in the outcome of the model 
without any further approximations other than the random phase approximation (RPA) together with the mean field 
approximation (MFA). 

With this in mind, our work improves, and at the same time confirms, the results of Ref. (36|, since we compute 
on the same footing both the meson masses and their effective chemical potentials, in the neutral phases. Moreover, 
at large enough density we find a new collective mode with K~ quantum number but with lower mass, standing just 
below the Landau damping threshold. This is an unexpected feature of the model, which comes out after the low 
density approximation is relaxed. Even if we cannot exclude the possibility that this new excitation it is just a model 
artifact, it should be stressed that it might modify the arguments about the critical baryon density for the onset 
of strangeness, and therefore it might even bring a drastic change in the conventional picture for the population of 
strangeness in the compact stars. 
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We also compare the NJL scenario studied here with that of a simple nuclear model in which kaon condensation 
is realized [HI, |3j| ; this comparison is useful since it suggests which are the possible directions to follow in order to 
reproduce, at least qualitatively, the scenario accepted by the nuclear matter physics community. For completeness 
we notice that other studies exist in which the kaon condensation scenario is not favored (40| . 

The plan of the paper is as follows: in section II we specify the lagrangian of the vector extended NJL model. 
In section III we derive the mean field effective action of the VENJL model as well as the effective action for the 
meson excitations at the second order in the meson fields. In section IV we show results for the restoration of the 
approximate chiral symmetry at finite baryon chemical potential. In section V we present a detailed discussion of the 
charged pscudoscalar modes in neutral and f3— equilibrated quark matter. In section VI, inspired by our results, we 
briefly discuss on the role that vector interaction might have on color superconductivity. Finally in section VII we 
summarize our results and draw our conclusions. 



II. THE VENJL MODEL 



In this work we are interested in neutral and /3-equilibrated quark matter at finite density. We work in the grand 
canonical ensemble formalism, introducing a chemical potential fi corresponding to the conserved baryon number. To 
be more specific we consider a system of u, d and s quarks at a finite chemical potential \i described by the lagrangian 

£ = [idfj,^ + ^7o] tl>f + Cmass 

f 

+ d + C 6 + C v + C A 

+ e [Wrf] e , (1) 

where ipf corresponds the quark field with flavor / (= 1,2,3 for u,d,s) and e denotes the electron field. We now 
specify each term in Eq. (flj). A sum over color indices is understood in Eq. (fTJ). 

Equilibrium under weak interactions d — > uev, s — > uev implies fi^ = fi 8 and /i^ = /i u + /i e (wc assume neu- 
trinos escape from matter therefore [i^ = 1 ). Moreover, in order to achieve electrical neutrality we introduce the 
Lagrange multiplier /iq associated to the total charge Q, add the term (J-qQ to the lagrangian |T]) and minimize the 
thermodynamic potential under the stationarity condition dfl/dfiq ~ 0. The total charge operator is given by 

2 11 

Q = -u f u - -tfd - -s f s - e f e , (2) 
o o o 

thus adding [iqQ to Eq. |T]) and recognizing that /xq = — fi e we get the chemical potentials of the quarks, namely 

M« = M - gMe , = Ms = M + g/^e • (3) 

Next we describe the other terms in our lagrangian ([1]). The mass term is 

C m ass = - ^mflpflpf (4) 
/ 

and m/ is the current mass. In this work we assume m u = m<j. The NJL four-fcrmion and six-fcrmion interaction 
Lagrangians are @, H, 

8 

£ 4 = [(fiXa^f + (ih^Y] , (5) 

a=Q 

C e = —K [detV>/(l +js)ipf> + det^/(l - 75)^/'] j ( 6 ) 
where A a are the Gell-Mann matrices in flavor space (Ao = \/2/3 1/) and the determinant is in flavor space as well. 



Some interesting possibilities may also arise when the neutrinos are trapped in the very early stage of the thermal evolution of proto- 
neutron stars before the deleptonization. The readers are reffered to l4lt l42l, |43|. |44| in this context. 
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Finally Cy + C-a denote the following U(3)y ® U(3)a invariant interaction term: 

8 

£v + C A = -G v Y,[(^^) 2 + (^j 5 Xa^) 2 ] , (7) 

where a summation over color and flavor is understood. In the above equation A a denote the same set of matrices in 
the flavor space introduced in Eq. JT]). Typical values of Gy in the vacuum are 0.2 < Gy jG < 3 fl6l. [l9l. l29l, l32l. |35| . 



III. EFFECTIVE POTENTIAL AND EFFECTIVE ACTION OF MESONS 



In this section we sketch the derivation of the mean field effective potential as well as of the effective action of the 
meson modes in the VENJL model. The derivation is done in some detail in order to easily compare our notations 
and results with the existing literature. We follow the conventions of Ref. Q- 

The easiest way to derive the meson propagators is the so called linear approach [5f| H|| Ht| ■ Firstly we write the 
t'Hooft term as an effective 4-fermion interaction: this is a well known procedure Q therefore we do not insist on its 
details. After this is achieved the scalar and pseudoscalar interaction is written as 
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c 



mixing 3 



(8) 



a=0 



where the term C r 



generates the mixing between it , rf and rf modes and is not important in this context since 



we focus on the charged meson modes. The effective coupling constants are defined as follows: 
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In the above equations Of = (//). Then we define the following meson fields: 

v£ = GvH 1 ^ > 
a£ = G v ^ 5 X a iP . 

The partition function of the model can be cast in the form [47| 



(9) 
(10) 

(11) 
(12) 

(13) 



(14) 
(15) 
(16) 
(17) 



VMVtpVtpe 



iS 



VM = VcffDirVvVa , 



with the bosonized action given by 



S 



d 4 x 4> [i^d - m + 2<fi a ip\ a ilj + 2ir a ^i"f b \ a ip - 2v^ipJnKip ~ ^t/jj^Xa^] ip 



d A x 



(-) 



K {+) 









Gy 



(18) 



(19) 



In Eq. I|19p we have not shown explicitly the contribution of the free electron gas (it will be inserted at the end of the 
calculation). At this stage the meson fields arc external fields (i.e. with no kinetic term). Their kinetic terms as well 
as their interactions will arise once the fermions are integrated out. Before performing this integral we notice that in 
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our application we expect condensation in some of the <p a and v£ channels. These condensations are related to chiral 
condensates of the three flavors and to fermion densities. Thus we define 



a = Va + $a 
v a 



a 



p 1 : 



(20) 
(21) 



with a = 0, ... ,8 and (S a ) = (p%) = 0. With this choice T>M = 'DSVifDp'Da and the physical meson fields will be 
identified with the fluctuations around their expectation value. In principle we should introduce expectation values for 
charged pseudoscalar mesons as well, since their in-medium energies can be lower than the threshold of condensation. 
However we firstly study the system without the assum ptio n of charged pseudoscalar meson condensation: as long 
as we exclude the possibility of first order transition [(381. l69l| , the evolution of their rest energies as a function of the 
baryon density will allow to establish wether condensation occurs or not (see next section for more details). 

The functional integration over VipVip in the partition function can be done exactly since the action (|19[) is quadratic 
in the fermion fields. Using textbook relations we get 



VMe iS[^w,v,a] ( 

with the effective action S[<fi, w, v, a] given by 

S[<j>, 7T, v, a] = -zTrlog [S^ 1 + 26 a X a + 2nr a j 5 X a - 2p^X a - 2a^^ 5 X a ] 
d A x 



(22) 



d 4 a 



V£ + pi + 2V aPa + ai 



< 1 


Ki- 




el 2 - 





and 



5*0 1 = ~ m + M7o + 2a a X a - 2V£-y li X a . 

For seek of compactness we introduce the total fluctuating meson field 

M a = S a + iir a ^ 5 - p^ - 0^7^75 . 
Expanding the log term we finally cast the effective action in the form 

S[<l>,ir,v,a] = Smf-iY,— Tr[(S 2X a M a y 



(23) 

(24) 
(25) 



n=l 



d 4 x 
d 4 x 



5 2 a + 25 a a a , n 2 a 



Pl + 2V aPa + d\ 



Gy 



where Smf denotes the contribution of the effective action that does not depend on the fluctuating fields, 



S 



AIF 



-iTrlogS'o 



d 4 x 



d 4 x 



V 



G\ 



(26) 



(27) 



The sum in Eq. (|26p is a loop expansion that generates the kinetic terms for the mesons as well as their interactions. 
We now examine Smf and this loop expansion separately. 



A. The mean field effective potential 

We now focus on the mean field term in Eq. (|27|) . Firstly we have to specify the expectation values in Eqs. ([20]) 
and (jnj). We assume that in the ground state condensation occurs only in the channels (//) and (/'/) (the- latter 
being relevant at finite quark density). Here / denotes the quark with flavor /. This implies that in Eqs. (|2U)) and (JHJ) 
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only the terms with a = 0, 3, 8 and p = survive, the reminders being zero. It is easy to show that S 1 takes the 
simple form 

/ K o o \ 

S^ 1 = h d (28) 

V o o /i s J 

where 

ft/ = (p + A*/ - 4Gyp/)7o - p • 7 - M/ , (29) 

and 

Mf=mf — iGcrf +2 K af+iaf+2 (30) 

denotes the mean field (or constituent) quark mass. In the above equation 

ct/ = -iN c txSf , (31) 

denotes the chiral condensate of the flavor /, where Sf is the propagator of the quark of flavor /, N c is the number 
of colors, and the trace is on spinor indices only. Also we have defined a± = a u , a 5 = ad- Analogously pf = (f f) 
denotes the number density of /. 

From Eq. (|29p the role of the vector interaction in the mean field approximation is clear: it shifts the quark chemical 
potentials to 

lt f ->Hf- 4G vPf . (32) 

Thus depending on the sign of Gy the chemical potentials will be shifted either upwards or downwards once the 
number density pf 7^ 0. 

We notice that the flavor channels Ao, A3 and As in the vector interaction generate coupling of the quark currents to 
the total quark number density p u + pd + p s , the isospin density p u — pd and to the hypcrcharge density p u + pd — 2p s 
respectively. This is similar to what happens in the rclativistic mean field nuclear models [38|, where the coupling 
of the isospin nuclear current to the isovector p M = (PijP^iPs) an< i of the nuclear current to the vector uj^ meson 
gives rise at the mean field level to expectation values of the zeroth component of the u* 1 and p 1 ^ fields, which result 
proportional respectively to the total baryon and to the isospin densities. We simply mention here that one cannot 
identify the vectors of the channels Ao, A3 and As with the physical u), po and <j> fields since a mixing among them 
occurs in the vacuum as well as at finite density. Keeping this in mind, with an abuse of notation we define 

= Gv(u^u + cPd + Js) = G v (Pu + Pd + Ps) , (33) 
(p°) = G v (u*u - dU) = G v (Pu - Pd) , (34) 
(0°) =G v {iJu + dU-2sU) =G v (p u + pd-2p s ) , (35) 



or equivalcntly 
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Pu = 


6Gy 




1 


Pd = 


6Gy 
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Ps = 


6Gy 



(2(^°) + 3( P3 >) + (0°)) , (36) 
(2( w °)-3(p§) + (0 )) , (37) 
(2<u,°) - 2(0°)) . (38) 



The above equations allow to rephrase the effective chemical potential in Eq. (|32|) in terms of the expectation values 
of the vector meson fields. 
Finally writing 

e-? vn = Z (39) 
with Z given in Eq. ([22]) we get the mean field effective potential: 

Qmf = -4Ka u <T d a s + 2G ^ a)- 2G V J2 Pf + V ^ > ( 40 ) 

f—u,d,s f—u,d,s 
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with 



-T 



+ 00 

E 



d 3 p 



Tr log 



h u 
/i d 
0ft, 



(41) 



In the above equations we have introduced a finite temperature T in order to easily handle infrared divergencies that 
arise when the chemical potential of the flavor / is larger than its mean field mass Mf. At the end of the calculations 
we make the limit T -> 0+. 

The ground state of the model at hand is characterized by the numerical values of the condensates and of the quark 
number densities. We determine the values of the chiral condensates by looking for the global minima of the total 
effective potential: 



f2 — Qmf — 



Me 



12tt 2 



(42) 



where the second addendum on the r.h.s is the thermodynamic potential of the free electron gas that will neutralize 
the net quark charge. At a fixed value of p we choose the value p e of p e which neutralizes the system, which is defined 
by the stationarity condition 







dfl e 



(43) 



The complete numerical strategy will be described in the next section. 



B. Effective action of meson fluctuations at the second order 



We now discuss the action of meson fluctuations in Eq. (|26p . The term with n = 1, when summed to the linear 
terms in the second and third line of the same equation, gives rise to terms of the form 8 a F a or p a G a , where F a and 
G a schematically denote the gap equations whose solutions determine the physical values of the chiral condensate and 
the quark number densities [471 ]. Since for these values both F a and G a vanish, the linear terms do not appear in the 
effective action. 

Next we consider the term with n = 2. This gives rise to the kinetic terms of the meson fields as well as to the 
mixing of some of the modes. Since the linear terms vanish we can write the effective action at the second order as 



S[S,n,p,a] = 2i d 4 x d^y Ti[So(x,y)X a M a (y)So(y,x)X b Mb(x)] 



d A x 



5a{x) 2 TT a {x) 



2 1 



K (+) 



+ / cf X 



p a (x) 2 + a a {xf 
Gv 



(44) 



where we have not written explicitly the terms proportional to a 2 , V 2 since they contribute only to Smf- In the 
above equation the trace is understood over color, flavor and Dirac indices. 

By means of the loop expansion in Eq. (|44|) we could study every meson fluctuation at the second order in the 
meson field about the mean field solution. However, in this paper we focus our attention on the charged pseudoscalar 
modes, reserving a more complete study to a next paper. 

The propagator of these modes is easily obtained from Eq. (j44|) . We are interested in the rest energies of the meson 
modes. We have verified by a direct calculation that there is no mixing between the pseudoscalar and the longitudinal 
component of the pseudovector mode in the case we put the spatial momentum of the meson Q = 0. Moreover parity 
conservation in strong interactions is enough to ensure that mixing at the second order in the fields does not arise 
between pseudoscalar mesons and the vector or the scalar ones. Thus the remaining term to be computed is that of 
order 0(Tr a iTb) of Eq. (|44|) . Fourier transforming both the fields and the propagators we have 



(?) 



( + ) 



2i 



_dH_ 



Tr {S (£h 5 X a S (£ + qh 5 X b ) 



(45) 



Since at the second order there is no mixing among the pseudoscalar and the other fields, at least in the rest system of 
the mesons, we identify the matrix in parenthesis in Eq. (|45[) with half of the inverse meson propagator in momentum 
space: 

I -UtTt [So(£)j 5 X a S (£ + qh 5 X b ] . (46) 



1 



§A4<?) _1 



Sab 



K. 



(+) 
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Consider for example the charged kaons whose fields are defined as 



K± = ^7? ' (47) 



then using D 4 l = D 5 ^ = D and D 4 § = — -D54 = iSD 1 it is easy to show that 

the condition D _1 (go = ^k,Q = 0) — <5 ~ 1 ( go = W A',<7 = 0) = with cj^ being positive defines the rest energy of 
charged K + . Since D~ 1 (—qo,q) = D^ 1 (qo,q) and 6D~ 1 (—qo,q) = —5D^ 1 (qo,q) can be easily shown, it follows that 
the negative energy solution to D^ 1 — SD^ 1 = corresponds to the positive energy solution to I? -1 + SD^ 1 = 
indicating that we need only one equation D _1 — SD -1 = to extract both the K~ and K + energies. Similar results 
hold for the remaining mesons. 

IV. THE NEUTRAL GROUND STATE OF THE VENJL MODEL 

In this section we discuss the neutral and /3-equilibratcd ground state of the VENJL model at zero temperature, 
obtained by the procedure of minimization of the effective potential Q under the neutrality condition (|43[) . We present 
results for the following set of parameters [48[ : 

m u = m d = 5.5 MeV , (48) 
m s = 140.7 MeV , (49) 
A = 602.3 MeV , (50) 
G = 1.835/A 2 , (51) 
K = 12.36/A 5 , (52) 

but the results are qualitatively similar to those that we obtain with the different set of Ref. [49[. Firstly we are 
interested in the effect of the vector interaction on the neutral ground state of the model. 

The value of Gy can (and should) be fixed in the vacuum in order to reproduce the spectrum of the vector mesons, 
see for example [l8j]. However, it is not clear how large is the modification of Gy in the medium [3l| . therefore instead 
to fix it by the meson spectrum in the vacuum we treat it as a free parameter in order to grasp its effects on the chiral 
restoration transition at finite density, as well as on the in-medium meson properties. We measure the strength of the 
vector interaction in terms of the scalar coupling G and introduce the ratio 

rv = % - (53) 

In the vacuum we find ry sa 0.6 in order to reproduce the neutral p meson mass m p = 775 MeV. We notice that because 
of the shift in Eq. (|29|) the quark number density for the flavor / in the mean field approximation is sclf-consistcntly 
defined by the equation 

PS = [0* " - IGvPf) 2 ~ M]] A/ 6 [((i - Qfi e - AGypff - M 2 f ] , (54) 

o7T 

where Mf denotes the in-medium quark mass. The r.h.s. of the above equation is nothing but the number density of 
a degenerate Fermi gas with a density dependent Fermi momentum given by 

[{v-Qtie-lGvpff-M}} 1 ' 2 . 

It is interesting to note that Eq. (|54[) is equivalent to the partial-derivative of the effective potential with respect to 
Pf which is related to condensations of u>, po and <j> through Eqs. (|33H35|) i.e., 

^ = 0. (55) 
dp f 
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In this case £Imf should be regarded as a function of eight variables, i.e., {cr/, pj, p e , p} [3(|. It turns out that the 
ground state is realized as the maximum with respect to variables pf when Gy > (repulsive vector interaction). 
However it does not immediately mean the ground state is unstable against the vector meson fluctuation. Actually 
the vector meson mode is not tachyonic as we always find a positive and real mass for p meson. This means the 
curvature of the effective potential cannot help us to judge if the system is unstable against small fluctuations in the 
vector channel. The total charge is given by 

2 11 

= -^Pu - ^Pd - -^Ps - Pe , (56) 

with quark number densities identified with the solutions of Eq. (|54"|) . 

Since the number densities in the case ry ^ have to be computed self-consistently by virtue of Eq. (|54|) the 
numerical calculations are slightly more involved than for the case ry = 0. Our numerical strategy is as follows: for 
each value of p we solve Eq. (|54[) for each flavor; this allows to define three numerical functions Pf = pf{p e ,Mf) 1 
one for each flavor /. Moreover we define a fourth numerical function which relates the physical electron chemical 
potential (namely the value of p e that corresponds to a vanishing total electric charge) to the in-medium quark masses. 
This is achieved by means of the total charge Eq. (|5"6")) and densities defined above. We insert the four functions into 
the effective potential, that now depends only on the in-medium masses (or equivalently on the chiral condensates). 
The global minimum of the effective potential in the space (M u , Md, M s ) gives the physical values of the in-medium 
masses. Once these are known we can step backward and compute the numerical value of p e and of the number 
densities. 

The phase diagram of the model can be described equivalently in terms either of the chiral condensates or of the 
mean field quark masses. We prefer the latter point of view since it is physically more intuitive. From the qualitative 
point of view the behavior of the constituent quark masses as a function of p is similar in the cases ry = and 
ry 7^ 0. In both cases there exist a critical value of p such that Mf for / =u,d suddenly decreases from its vacuum 
value Mq « 367 MeV to a lower one, typically a few tens of MeV. However, the order of the transition depends on 
the value of ry. For ry < we find a first order transition as in the case ry = 0. In this case the critical chemical 
potential is lowered with respect to the case ry = 0. The first order character of the transition remains even for small 
and positive values of ry. At ry = r v ss 0.5 the first order transition becomes a crossover. In this case we identify 
the crossover point with the value of p where \dM u /dp\ is maximum. The magnitude of \dM u /dp\ at the critical 
point is lowered as ry is increased above the ry. indicating that the crossover becomes more smooth. These results 
are in qualitative agreement with those obtained within the non neutral and non /3-equilibrated VENJL [hJ, [2§| ; they 
are also in agreement with the results obtained within the three flavor PNJL model in Ref. [3l|, where it is shown 
that the first order line disappears from the T — p phase diagram when Gy is increased above a critical value. We 
summarize these results in Fig. [1] where we plot M u and M s against p for three representative values of ry (left 
panel), as well as the critical chemical potential as a function of ry (right panel). For comparison we plot, by dashed 
line in the right panel of Fig. [TJ the critical chemical potential computed without requiring the neutrality condition 
and by setting p e = 0. We notice that the neutrality condition has an effect similar to that of the repulsive vector 
interaction: it stabilizes the chiral symmetry broken phase with respect to the restored phase, as it increases the 
critical chemical potential upwards. The same effect can also be induced by increasing the strength of repulsive vector 
interaction. Interestingly enough, the value of critical chemical potential at which the first order phase transition 
turns into a smooth crossover is not so much affected by the neutrality condition. However, the required value of the 
strength of repulsive vector interaction is significantly reduced if the condition of neutrality is taken into account. As 
a consequence, the neutrality constraint actually helps to make the transition smoother at fixed vector coupling. 

In Fig. [5] we plot the baryon density pb defined as 

Pb = r i iPu + Pd + Ps) , (57) 

as a function of p for four different values of the ratio ry = Gy/G. The densities of each flavor are computed by 
means of Eq. (|54[) with the values of p e and Mf obtained by the minimization procedure. In the figure the dot-dashed 
line corresponds to ry = —0.2, dashed line to ry = 0, solid line to ry = 1.1. An interesting consequence of the change 
from first order phase transition to crossover into the approximate chiral restored phase when ry > 0.5 is that the 
baryon density as a function of p, is a continuous function of p. Thus the system smoothly passes from a dilute Fermi 
gas to a dense one, the size of the smoothness depending on the precise value of ry. This does not happen when the 
transition is of first order. 

In the left panel of Fig. [3] we plot the electron chemical potential of neutral quark matter as a function of p for 
some value of ry. In the right panel wc plot the same chemical potential against the baryon density ps- The latter 
plot is obtained by assembling data from the left panel and from Fig. O It is interesting to notice that at a given 
value of p, the larger the magnitude of ry the larger the numerical value of p e . For example at p = 440 MeV wc 
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FIG. 1: Left panel: Mean field up quark mass (thick lines) and strange quark mass (thin lines) against p for three different 
values of the ratio ry = Gy/G in neutral and /3-equilibrated quark matter: dot-dashed line corresponds to ry = —0.2, dashed 
line to rv = 0, solid line to ry — +1.1. Right panel: critical value of p for restoration of the approximate chiral symmetry as a 
function of rv. Bold line denotes first order phase transition, solid thin line corresponds to a smooth crossover. For comparison, 
the same quantity computed without requiring electrical neutrality is shown by dashed line. 




FIG. 2: Baryon density in neutral quark matter, in units of the saturation density po = 0.16 fin 3 , against p for four different 
values of the ratio ry = Gy/G: dot-dashed line corresponds to ry — —0.2, dashed line to ry — and solid line to ry — 1.1. 



find p e (ry — 0) ~ 95 MeV, to be compared with /i e (ry = +1.1) w 120 MeV. For comparison, at the same value of 
p we find p e (ry = —0.2) sa 80 MeV. The change of the sign of the slope of the curves at large p (or ps) occurs in 
correspondence of the condition p s > (at lower values of p we find p s =0), which implies that strange quarks take 
a role in the neutralization of the system and a less number of electrons is needed. 



V. MESON ENERGIES AND (ABSENCE OF) CONDENSATION 



In this section we compute the pscudoscalar meson energies as a function of the mean quark chemical potential p 
at T = 0. To achieve this result we solve the pole equation in the rest frame for the appropriate channel as discussed 
in a previous section. We focus on the charged modes here because they are interesting in the context of meson 
condensation driven by the electron chemical potential [H, [H, [H, [13, El, [69|, [lO, [7l|, [ll, E3] • 
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FIG. 3: Left panel: electron chemical potential in neutral quark matter against fj, for four different values of the ratio rv = 
Gv/G: dot-dashed line corresponds to rv = —0.2, dashed line to rv = 0, solid line to rv = 1.1. Right panel: electron chemical 
potential against baryon density. The line sketch is as in the left panel. The change of the sign of the slope of the curves at 
large fi (or p_e) occurs in correspondence of the creation of strange quark Fermi spheres (/i s > M 3 ). 



A. Meson energies and in-medium meson propagator 

Keeping into account the results discussed in the previous section (see Eq. (|45jl ) we write the equation for the 
energy of the charged kaons as 



ftF K ±{uj) = 0, 
F K± (u) ee l-2ivi +) n K± (^,0), 



where the polarization function is defined as 

n K ±(u;,Q) = -2iN c J -^tv 



05 



r75 



h u (p) h s (p + Q) 



(58) 



(59) 



where Q = (lo,Q) with Q being the three momentum of Kaon, and denotes the real part. In the above equation 
the trace is understood on Dirac indices only. Performing the trace the pole equation (|58|) in the rest frame reads: 



243? 



d 3 p 



T E 



n— — 00 



+ p 2 - (iw n + Hu)(iun + i H~ Ms) 



1 



[(iw„ + [iu? P 2 M 2 } [{iuj n + in m + f i s ) 2 -p 2 ~ M 2 } 2K ( + ] 



(60) 



Once again we have introduced a finite temperature T in order to handle infrared divergencies that arise when the 
chemical potential of the flavor / equals its mean field mass Mf. At the end of the calculation we put T — > + . In 
the above equation f2 m = nTn is the boson Matsubara frequency. 

The retarded real time propagator is obtained via iQ m — ► cu + i0 + after summation over fermion Matsubara 
frequencies. The result turns out to depend on the external energy u> only on the combination 



) = w + [l s - n u = UJ + n e + 4:G v (p v 



(61) 



Thus also the kaon lagrangian in momentum space depends only on Qq. This implies that in the derivative expansion 
one can build only terms that contain the covariant derivative 



iD Q EE id + He+ 4Gv(pu ~Ps) ■ (62) 

We make analytically the sum over Matsubara frequencies in Eq. (foT))) . After continuation to real external energies 
we take the limit T — > + , and we get an expression that depends on energy only through Q . We do analytically the 
integral over \p\. Finally we solve Eq. ([BU)) in the variable Qq. The solution of the pole equation in the Qq variable 
defines the in-medium kaon mass m* K± . Within our convention in Eq. (|59|) the positive (respectively negative) energy 
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solution of Eq. ([60]) . Qq = m* K _ (respectively Qo = —m* K+ ) correspond to the K (respectively K + ) in-medium 
mass. Once the pole equation is solved at a given value of p, we obtain the in- medium energies trivially from Eq. (|6ip : 

= Me + 4G v (p u - p a ) +m* K+ , (63) 
Uk- = -A* e - ^G v (p u - p a ) + m* K - . (64) 

At Gy — the chemical potential felt by kaons is — qp e where q denotes the electric charge (in units of the electron 
charge) of the meson. When Gy is switched on, the kaon chemical potential is shifted from — qp e by virtue of the 
strong interactions. As a matter of fact, since the quark densities are related to the expectation values of the p, lo and 
<\> mesons, see Eqs. (|36|) and ([38]). the Gy terms in the meson energies can be interpreted as due to the interactions 
of kaons with the aforementioned vector mesons (there is also a modification of m* K when Gy is switched on). We 
stress however that these interactions do not exist at the tree level within the model we study in this work, but arise 
only as loop effects: the quarks interact with the expectation values of the vectors and the quark loops, that generate 
the mass and the kinetic terms of the pseudoscalar mesons, will depend on these expectation values. Another source 
of interactions among pseudoscalar and vectors (as well as axial-vectors) are the terms of the cubic order in meson 
fields in the loop expansion of Eq. (|26p . that we do not consider in this study for simplicity and will be the subject 
of a next paper. 

The other pseudoscalar channels are treated in a similar way. The pole equation for the charged pions reads 

24S{? f d 3 p MyMd + P 2 - (iWn + Pu)(iUn + i^m + Pd) _ 1 

J (2tt)3 [{iu> n + p u f p 2 - Ml] [(iw n + ifi ro + pd? ~P 2 Ml] ~ 2K[ +) ' 1 ' 

and in this case the polarization tensor depends on external energy only through the combination 

Qo = u + p d - Pu = OJ + p e + 4Gy {p u - p d ) . (66) 



In this way we have 



^7r+ = Me + 4Gy(/>„ - p d ) + m*+ , (67) 
u n - = -p e -4G v (p u - p d )+m*- . (68) 



Similarly, the equation for the neutral kaons is given by 
d 3 p ^ M d M s +p 2 - 

( 2 ^) 3 nhoo ^ + -P 2 - M l\ [fan + lSl m + p s f -pi- M|] 2l4 +) 



24 sp f d 3 p M d M s + p 2 - {iuj n + p d )(iun + i^m + Ps) = 1 (( , q) 



with 



Q = uj + p s -p d = u} + AG v {pd - Ps) ■ (70) 

The in-medium energies are given by 

u> K o = AG v (p d - p s ) +m* K0 , (71) 

u k° = -4Gy(p<j- Ps) + m* RQ . (72) 

We solve Eqs. ([63)) and (f6"5)) by using the values of <jf and p e that correspond to the neutral global minimum 

of the effective potential (|4"2"|) . In Fig. 0] we plot the solutions of the real part of the pole equations for charged kaons 
(upper panels), corresponding to Eqs. (|63|) . (f64|k neutral kaons (middle panels), corresponding to Eqs. ([67]) , ([68]) : and 
charged pions (lower panels), corresponding to Eqs. |[7T]) . (|72p . They are plotted against p for ry = (left panels) 
and ry = +1.1 (right panels). Results for intermediate values of ry do not differ qualitatively from those shown in 
the figure. Also, for negative values of ry the difference with the case ry = is that the location of the discontinuity 
in the in-medium energies is shifted to lower values of the quark chemical potential. 

The plots of Fig. 0] are interesting for several reasons. A general feature of the vector interaction is the stabilization 
of the charged pseudoscalar states: the comparison at the same chemical potential of the in-medium energies of a 
meson with and without keeping into account the repulsive vector interaction shows that they are higher in the latter 
case than in the former one. Secondly, we notice that the effect of the vector interaction on the K~ and K° in-medium 
energies is to lower them in the p regime where chiral symmetry is approximately restored. For example at p = 460 
MeV we find uj k - ps 406 MeV for ry = 0, to be compared with uo K - sa 290 MeV for ry = 1. This is in part due to 
a mild lowering of m K as Gy is increased, but mainly to the effective chemical potential arising from interaction of 
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FIG. 4: Solutions of the pole equations in neutral quark matter for charged kaons (upper panel), neutral kaons (middle panel) 
and charged pions (lower panel) as a function of for two different values of ry. 



K~ and K° with the expectation values of the vector meson fields via the quasi-quarks in the loop, as is clear in the 
expressions, Eq. and Eq. (f?2")) . In the case of K~ there is still a further lowering of the energy due to the larger 
value of /Lt e , see Fig. [3] at /i = 460 McV we find /i e r* 95 McV for ry = 0, to be compared with [i e w 120 MeV for 
TV = +1.1. 

For completeness, in Fig. [5] we plot the separate contributions to the in-medium Kaon energies. The solid lines 
indicated by ttlk*_ express the in-medium Kaon masses while /_t e is electron chemical potential. The in-medium kaon 
mass ra if* should be compared with fj, s — fj, u = fj, e + 4Gy( J o u — p s ), which is effective chemical potential felt by K~: 
it is shifted from the bare charge chemical potential /j, e means of the vector interaction. Both the electron chemical 
potential and the vector interaction tend to lower the in-medium K~ energy, even though the lowering is not enough 
for the kaon condensate to form. 

The shape of the inverse K~ propagator at rest as a function of energy and /x has an interesting feature. In Fig. [6] 
we plot the real and the imaginary part of F K ± (w) = 1 — 2K^H K ± (ui, q = 0) for ry = +1.1 and two values of fi. They 
are represented respectively by the solid and the dashed line. At /i = we have 'SiF K ±(niK) — 31Fk± (— ttik) = 
and QF K ±(u) = for \w\ < M u + M s . Here ttik denotes the vacuum kaon mass. As fi is increased above the chiral 
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FIG. 5: The various contributions to in-medium Kaon energies. Effective in-medium Kaon masses rriK* ± , the electric chemical 
potential [i e and /j, a — [i u as a function of fi. mK ± * ± (fJ-s — Hu) correspond to the two curves in the top-right figure of Fig. [4] 
The deviation of fj, s — fi u from fj, e comes from the repulsive vector interaction. 




FIG. 6: Upper panel: Real and (muinus) imaginary part of F K ± (lu) si — 2K^Tl K ± (lj, q = 0) as a function of u for Gv = 
+1.1(7, for two different values of /j,. Solid line corresponds to the real part, dashed line to the imaginary part. Lower 
panel: Corresponding spectral function of Kaonic propagator p(u)) = Q (1/F k ±(lu)). 



transition then Fk± (w) develops two singularities at intermediate values of u>, see the right panel of Fig. [6l Moreover 
two new zeros appear. Thus in principle at large fi four kaon modes appear. The solutions we have chosen to draw 
Fig. [S]are those with the larger magnitude. 

In order to make clear the physical content of the intermediate singularities at u> = u)\, L02, we will take a detailed 
look at the imaginary part of F K ± (lu) from now on. It is possible to find some analytic formula for the imaginary part 
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FIG. 7: The energy of the collective modes in the (w-fJ.) plane. The new collective mode with quantum number K~ induced 
by the Landau damping. The shaded area corresponds to the region where the responce function suffers from imaginary part 
due to the Landau damping. u) K - and lo k + are the same as two curves in the top-right figure of Fig. [4] The points where the 
solid lines turn into the dashed ones indicate the thresholds for the continuum: across the points, the modes become no longer 
stable decaying into the constituents. 



which is physically related to the K decay rate in the medium. By use of the formula $s(l/ (u}+x+i8)) = —iir5(uj + x), 
we find 
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(a) 
(b) 
(c) 
(d) 



(73) 



where Qo is again the shifted energy denned by Eq. (fBTj) , E u = M£ + p 2 and E s = \J M 2 + p 2 are the u and s 
quark energies, and f's are the Fermi blocking factors defined by /j = 6(fXi — Ei) for the quark fermi distribution, 
and fa = 9(—fii — Ej) for the antiquark Fermi distribution. Since the momentum integral is definitely convergent, we 
have removed the cutoff from the evaluation of the imaginary part. As easily guessed from the blocking factors and 
delta function guaranteeing the energy conservation, (a) expresses the two body decay from K~ to (su), (6) does the 
Landau damping process of K~ absorbing onshcll u quark decaying into s quark, (c) represents the inverse Landau 
damping of K~ absorbing anti s quark decaying into anti u quark, and (d) is for the decay of K~ by absorbing u 
quark and anti s quark. The p-integral is extremely trivial due to delta functions, yielding 



%F K ±(tu) 
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(a) 
(b) 

(c) 
(d) 



(74) 



where E 



<*) - 



Mf + an( i {P*?P*} are now to be determined through the condition of energy conservation, 



E* + E* = \Qo\ for |Q | > M u + M s 
Et-E* = [Q \ for |Q | <M S -M U 



(75) 
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each of which has a unique solution at fixed Qq in its domain of definition. It turns out both have the same functional 
dependence on Q 0l i.e., 



P ^ = m ■ (76) 

We now consider in which condition the Landau damping process labeled by (b): u + K~ — > s, (and the inverse 
production process (s — ► it + A' - )) is possible. As demonstrated above, the energy conservation £7* = Q + E* has a 
unique solution when Qq is in the interval: 

< Q < M s - M u . (77) 

From Eq. ((75]), we see is a decreasing function of Qo, and when Qo approaches zero diverges. In addition to the 
energy conservation, the Pauli principle should be also satisfied: when p+ > \/ — M 2 , or cquivallcntly E* > p u , 
there is no u quark avaiable to contribute to the absorption process u + K~ — > s. This can be also seen in the 
blocking factor /„ — f a in (b) of Eq. (|T4[) : This factor comes from the sum of two blocking factors, f u (l — f s ) for the 
decay, and — / s (l — /„) for the corresponding creation process, s — > u + K~. Then in order for the total decay rate 
to be nonzero, either {f u — 1, f s — 0} or {f u = 0, f a = 1} should hold. In our case M s > M u is always realized so 
we need to consider only the former condition, that is, pf s < P* < PFu, where the Fermi momentum is defined as 
PFi = v max.(0, pf — M 2 ). Solving this condition together under the kinematic constraint Eq. (|77[) . wc obtain the 
condition: y/Mf + p 2 Fu - yj M% + p 2 Fu < Q < \J M 2 + p 2 Fs - y/M% + p Fs , which, combined with Q = u + p, s - fi u , 
translates into 



w l = V M s + Pfu _ A*s _ max.(M u - /z„, 0), 

(78) 

^2 = ~\JMI + pp s +(i u + max.(M 5 - fi s , 0). 



We have checked that these equations reasonably reproduce the numerical values of {wi,^} in the right panel of 
Fig. [5] We conclude that the imaginary part in the range lo\ < lu < u>2 is non-vanishing due to the Landau damping 
of K~ by absorption of on-shcll u quark in the system. What is more intriguing is that below this Landau damping 
threshold, there appears a new pole implying a new collective mode with the same quantum number K~ as shown 
in Fig. [6j We have checked that this pole appears only at high density, fi < 420 McV, as depicted in Fig. [Jj At 
the moment we cannot discern if this new light mode is an artifact of the NJL model itself, as well as of the RPA 
approximation. However, if this is not the case, it is of great interest since it might change the low temperature 
behavior of the thermodynamics (e.g., strangeness population) in a compact star. We leave the detailed investigation 
of the nature of this new pole to a forthcoming paper. 

A third point that deserves discussion is the missed symmetry in the K° /K° spectrum. This splitting is present 
even in the case Gy = 0. It is due in this case to the fact that K° ~ (ds) and A'o ~ (sd): onshcll d quarks in the 
system tend to reduce the binding energy for Kq state by Pauli blocking. This is the reason why we see that Mk is 
higher than Kg . The splitting is more pronounced by the vector interaction, see right panel in Fig. [4] because of the 
induced effective chemical potential, Eqs. (fTTj) and ([72"]) . 

We finally comment on the absence of kaon condensation in our model. It is well established that at in the range 
of baryon densities (2.3 — 5)po and for neutral and /3-equilibrated nuclear matter, kaon condensation occurs, see for 
example [50l . l5ll . |52| and references therein. The main results of (50l . l5ll . |52| are in agreement with those that can 
be obtained within a simple relativistic mean field model [38l |39| . In all of these models the numerical value of /i e 
at a given baryon density ps is larger (at least a factor of 2) than the /x e we obtain in our calculations. Also the 
effective potential felt by kaons and due to interactions with p and lu mesons is enhanced in neutral nuclear matter, 
while in the model we study here it has only a mild dependence on pb- We guess that these two factors are the 
main source of difference between our results and the nuclear matter literature. At the moment we have not yet 
found a solution for this problem. A first possibility could be the hadronization of the NJL model [2l[, that would 
allow to describe the system in terms of protons and neutrons instead of quarks. Probably the neutralization and 
/3-equilibration of this system would result in a larger value of p e . Secondly, the expansion of the effective action of 
the meson fields to higher orders would introduce in the theory the direct coupling of kaons with lu and p mesons. 
However a dimensional analysis shows that these contributions are paramctrically suppressed by powers of (luq) /tuk, 
thus they should not be important. In order to verify this a direct calculation of the relevant diagrams should be 
performed. Thirdly, we cannot exclude that multi-quark interactions, that are known to be useful in the stabilization 
of nuclear matter 0, [3^, EH, can play a role in increasing the strength of the effective kaon potentials. Tensor 
interactions might also lead to an enhancement of the attractive interaction felt by the kaons in the medium. Last 
(but not the least) there is the possibility to include a density dependence of Gy via a fit to experimental data on 
the lu meson mass [33| . Some of these aspects are under current investigation. 
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FIG. 8: Quark number densities against p for ry = (left panel) and ry = +1.1 (right panel). 

VI. VECTOR INTERACTION AND COLOR SUPERCONDUCTIVITY 

Wc briefly comment on the effect that our results might have on color superconductivity [Hi], [H, HH, [fi3|. At 
large baryon density color superconductive phases, neglected in this paper for simplicity, could play a relevant role 
in the determination of the ground state. In the two flavor case the interactions in the quark-antiquark and diquark 
channels have been considered for example in Ref [35| . In order to understand how vector interaction can interplay 
with superconductivity in the three flavor model under investigation here, we compare the various quark densities at 
H = 550 MeV for r v = and r v = 1.1, see Fig. M 

In the case ry = we find p u ~ 12.9po, pd ~ 14.2/50, Ps ~ 11.6po; the in-medium quark masses ar M u w M ( i w 7 
MeV and M s « 197 MeV; the Fermi momenta of the different flavors arc p% w 538 MeVniw 556 MeV and p s F w 520 
MeV. It is well established that in these condition quarks condense to the CFL state |48l 158L . On the other hand 
for the case ry — 1.1 we find p u ~ 4.35po 7 pd ~ 5.95po, Ps ~ 2.71po; the in-medium quark masses ar M u m Md ~ 20 
MeV and M s « 390 MeV; the Fermi momenta of the different flavors are p F ps 374 MeV, ps 415 MeV and p s F ps 320 
MeV. It is not clear in this case which superconductive state could be favored, if any, since a self consistent calculation 
in color superconductivity that keeps into account the vector interaction is missing in the three flavor case. 

If wc suppose that color superconductivity sets in, in order to guess which state could be a good candidate in these 
conditions wc compute the ratio M%/fi. We find M%/(i ps 276 MeV. In these conditions, since p F — p F ps p F — p F « 
2(Pp- — Pr ) ~ 50 MeV, probably the three flavor crystalline LOFF state is the best candidate in the weak coupling 
regime [43, 0, . If the diquark coupling is high enough at this value of the chemical potential then a 2SC phase 
could be the ground state. Beside this there is probably room for spin one condensates. We also guess that because 
of the smooth crossover (see right panel of Fig. [8} and of the non trivial cou pling b etween the diquark and meson 
excitations a coexistence region might be created, as already discussed in Ref. [35|, [63|]. Of course these are just some 
of the possibilities, that should be either supported or not by a dynamical calculation. We finally notice that the our 
reasonings are not beyond the capabilities of the model even if we take p = 550 MeV since the Fermi momenta in the 
case ry — 1.1 are quite small compared to the ultraviolet cutoff A (they are smaller than the vacuum kaon mass). 

VII. CONCLUSIONS 

We have explored the consequence of the vector interaction in /3-equilibrated and neutral three flavor quark matter 
at finite density. Neutrality and /3-equilibrium are required to reproduce the conditions that coulb be realized in cold 
neutron stars. The extension of the NJL model to the VENJL one has not only an academic interest. As a matter of 
fact, in view of vector manifestation the vector mesons can play a relevant role in the restoration of chiral symmetry 
at large baryon density [HlHEl- Moreover, if we wish to describe the intermediate baryon density region of the NJL 
phase diagram in terms of a bosonized (and eventually hadronized) action, we need to take into account the effects 
of cj and p mesons exchange to be more realistic. 

Even if it is possible to choose the value of the vector coupling Gy in the vacuum, it is not clear what is its value in 
the medium. For this reason we have fixed the scalar coupling G and treated the ratio ry = Gy /G as a free parameter. 
We have found the interesting result that there exists a critical value of ry w 0.5 above which the approximate chiral 
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restoration becomes a crossover (at ry = the N JL model predicts a first order transition) . This result is summarized 
in Fig. [T] Our result is in agreement with Klimt, Lutz and Weise [l9[ and with Fukushima [111, where neutrality and 
/?-equilibrium were not required. 

We have also started a systematic study of the meson energies as a function of the baryon density and on the 
influence of the vector coupling, focusing in this explorative work on the pseudoscalar channels and leaving a more 
complete study to a next project. In particular our findings on the K~ energy show that kaon condensation is quit e 
hard to be realized within the present version of the model. This is disagreement with other results [H, H^, l5ll.l52j] , 
but we argue that the source of the disagreement is mainly the still poor description of the meson-quark interactions 
in matter at densities of the order of few times po that we have within the VENJL model. We have itemized some 
improvements that might lead to a better description of matter in the intermediate density region, and some of the 
aforementioned research lines are currently under investigation. 
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